Flat-Band Ferromagnetism in Organic Polymers Designed by a Computer Simulation 
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By coupling a first-principles, spin-density functional calculation with an exact diagonalization 
study of the Hubbard model, we have searched over various functional groups for the best case for 
the flat-band ferromagnetism proposed by R. Arita et al. [Phys. Rev. Lett. 88, 127202 (2002)] in 
organic polymers of five-membered rings. The original proposal (poly-aminotriazole) has turned out 
to be the best case among the materials examined, where the reason why this is so is identified here. 
We have also found that the ferromagnetism, originally proposed for the half-filled flat band, is stable 
even when the band filling is varied away from the half-filling. All these make the ferromagnetism 
proposed here more experimentally inviting. 

PACS numbers: 75.10.Lp, 71.20.Rv, 71.10Fd 



I. INTRODUCTION 

Materials design to realize desired properties in 
condensed-matter physics is becoming increasingly re- 
alistic. One important factor is, given high levels of 
computer performance, the computer simulation is an 
essential part in designing. Namely, computer simula- 
tions with first principles calculations should be impera- 
tive in narrowing down, or even pin-point, from a wide 
variety of candidate or designed materials for the inter- 
esting properties and functions we search for. Such an 
approach is especially promising for molecules and nano- 
structures. While the atom manipulation with scanning 
tunneling microscopy (STM ) 1 : 2 : 3 ! 4 provides new possibil- 
ities for materials design, the design of molecules, espe- 
cially polymers as we focus on here, should be fundamen- 
tal. 

In this context organic molecules and polymers are 
of special interest, since they have versatile structures 
and chemical properties that can be wider than in in- 
organic materials. Indeed, the discovery of a conduct- 
ing organic polymer by Shirakawa££ kicked off inten- 
sive studies that are paving a new way to adopt or- 
ganic polymers or oligomers in realizing various func- 
tions as molecular-electronics devices such as ficld-cffcct 
transistors^ or light-emitting diodes^. 

Organic ferromagnets have also attracted much at- 
tention as a challenging targels^^. In particular, or- 
ganic magnets consisting entirely of non-magnetic ele- 
ments is of fundamental as well as practical interests. 
Ordinary ferromagnets consisting of magnetic elements 
exploit electrons in d- or f-orbitals that are strongly inter- 
acting. Can magnetism arise in p-electron systems that 
are weakly interacting? One theoretical possibility is to 
apply the flat-band ferromagnetism proposed by Mielke 
and by Tasakii±*±£. This arises as an effect of electron- 
electron repulsion when the (one-electron) band structure 
contains a dispersionless band. The mechanism is inter- 
esting in many ways, but essential features are, first, the 
system is totally distinct from the "narrow-band limit" in 



textbooks, since the magnetism occurs when the transfer 
between different sites is finite. 

Second, this is a band ferromagnetism rather than a 
magnetism arising from the exchange interaction between 
localized spins. So, this mechanism should be a good can- 
didate for designing organic ferromagnets. While there 
are also other possible mechanisms such as the intra- 
molecular Hund's coupling for spins in p-orbitals&, the 
fiat band ferromagnetism is advantageous in that spins 
are carried by itinerant electrons which can be utilized 
for spin injectors in spintronics. 

Here a theoretical remark on the flat-band ferromag- 
netism is due. The flat-band ferromagnetism was first 
proposed by Liebji^ who considers bipartite lattices (con- 
sisting of two sublattices) that have different numbers of 
A sublattice sites Na and B sublattice sites Nb- Lieb 
proved that, when we switch on the electron-electron re- 
pulsion (assumed to be short-ranged, so that we take 
the Hubbard model), the ground state should then be 
ferromagnetic with the magnetization oc Na — Nb- Quan- 
tum chemically, this model contains Na — Nb non- 
bonding molecular orbitals, so is similar to Mataga's 
model. Shima and Aokiia then proposed a systematic 
way to realize such systems as superhoneycomb struc- 
tures. MielkeiSiifiiil and independently TasakiiS then 
constructed other flat-band ferromagnetism, which is dis- 
tinct from Lieb's in that the flat band is constructed from 
quantum mechanical interference between the nearest- 
neighbor and further transfers (so the system is neces- 
sarily non-bipartite) and the ground state is now ferro- 
magnetic. 

To be more precise, the lattice is required to sat- 
isfy what is called the "local connectivity condition'' . 
Namely, the flat band is by no means a sufficient con- 
dition for ferromagnetism, and the magnetism arises for 
special lattices on which adjacent "Wannier" orbitals 
have to overlap with each other no matter how they are 
combined to minimize each orbit size. So, while ordinar- 
ily a flat band implies a disjointed set of orbits, they are 
connected in Mielke- Tasaki lattices. This is intuitively 



2 



why spins align to lower the repulsion energy from Pauli's 
principle. Rigorous proof for the ferromagnetic ground 
state has been given for the half-filled flat band. 

In actual materials, it is very hard to realize perfectly 
flat bands, and to make them perfectly half-filled. For 
the former, a rigorous proof (for Mielke-Tasaki lattice)^ 
and numerical simulations*^! showed that a slight band 
dispersion does not destroy the magnetism. For the lat- 
ter, there are also a proo£ii and a numerical simulation** 
which revealed that small deviations from half-filling are 
allowable for the ferromagnetism. 

In the context of this background, what kind of or- 
ganic materials are promising for realization of the flat- 
band ferromagnetism? The hardest part is to realize 
the connectivity condition. In general we have to con- 
sider complex lattices such as Kagome, or lattices hav- 
ing distant-neighbor transfers. However, the condition 
is easier to satisfy on one-dimensional polymers. For 
example, Mielke-Tasaki model is realized as a chain of 
triangles i22i22*2^ So we can concentrate on polymers. 

Next comes the choice of the monomers that should be 
polymerized. Even-membered rings such as benzene are 
disadvantageous in that antiferromagnetic order tends to 
occur, which should dominate over the ferromagnetism. 
So we should opt for odd-membered rings, and the above 
example of the chain of triangles is indeed an example 
of this. Since triangular molecules are scarce, we can 
conclude that we should focus on polymers based on five- 
membered rings, such as polypyrrole, polythiophen, etc. 

For the computer design, the purpose of the present 
paper, there are two tasks: (i) to search for the materials 
that realize the flat band a la Mielke-Tasaki. Since the 
chain of odd-membered rings is by no means a sufficient 
(nor necessary) condition for the flat band, this is an 
important part. 

(ii) Next, we should do a first-principles (spin density 
functional) band calculations to confirm (a) the flatness 
of the band and the ferromagnetism in the ground state, 
and (b) whether the magnetism can indeed be interpreted 
as the Mielke-Tasaki mechanism, for which we have to 
evoke an electron-correlation calculation (exact diagonal- 
ization of the Hubbard model here), (c) We also address 
the question of how the departure from half-filling affects 
the magnetism in the flat band. 

For the first-principles calculation, we have employed 
the generalized gradient approximation 24 (GGA). De- 
pending on the purpose of each calculation, we have 
performed either the spin-unpolarized case based on 
the density functional theory (GGA-DFT) or the spin- 
polarized one based on the spin density functional the- 
ory (GGA-SDFT). We used the plane- wave-based ultra- 
soft pseudopotentials^Si. The energy cutoff was taken 
as 20.25 Rydberg. The convergence criterion of the ge- 
ometry optimization was that all of the forces acting on 
each atom were within 1 x 1CP 3 Hartree/a.u. 

The plan of this paper is as follows. In section [HI we 
describe the search for various materials and show some 
examples which exhibit possibilities of ferromagnetism. 




FIG. 1: Schematic view of (a) the eigenfunction that satisfies 
the connectivity condition and (b) the overlap of the adjacent 
eigenfunctions. 

In section HTT1 we show the result for the most promising 
material, poly-aminotriazole, and discuss how the theory 
for the flat-band ferromagnetism works on this particular 
material. In section lTvl we discuss how the ferromagnetic 
state is robust against the deviation of the band filling 
from the half-filling. 



II. EXAMINED MATERIALS 

A. Flat bands in the polymer of five-membered 
rings 

Let us first examine the simple tight-binding model, 

5 

H = -t ^ c \ c j + £ o X!" 4 + £lUl ' 

(id) *=2 

for a chain of five-membered rings, where we first assume 
that all the transfer integrals within the ring and those 
connecting rings have the same t, all the on-site energies 
are the same (eq = 0), except for the one at the top of 
the ring (ei). We first consider the half- filled case where 
one electron per site on average. This is a reasonable 
assumption when we consider 7r-orbital networks on such 
polymers. 

Let us start from an observation that Mielke-Tasaki's 
condition is satisfied when e\ = t (> 0), a not unrealistic 
condition. In this case, the third band with E = 0, i.e., 
the half-filled band, has no dispersion. This provides a 
heuristic example for the Mielke-Tasaki eigenfunctions: 
one can construct the most compact eigenfunction which 
consists of two rings as depicted in Fig. IHa) , where the 
amplitudes of the numbered sites are given as (1, 1, 0, 
-I, 0, -I, 0, 1, 0, -1). The fact that the amplitudes at the 
sites 5 and 7 are zero ensures the localized nature of this 
eigenfunction. Note that we are not displaying a part of 
a Bloch function, but the whole eigenfunction. One can 
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FIG. 2: (a) An atomic structure of polypyrrole. The structure 
extends along the right and left directions periodically, (b) 
The (GGA-DFT) band structure of polypyrrole. 



also construct another, linearly-independent eigenfunc- 
tion by simply shifting the wave function from one ring 
to the next. The eigenvalue of the new function is the 
same (E = 0). So the set of all these eigenfunctions can 
be a basis for the flat band. Because the eigenfunction is 
not confined in a single ring (unit cell) but extends over 
the two rings, two neighboring bases (eigenfunctions) al- 
ways overlap with each other as shown in Fig. ^b) . One 
cannot remove such overlaps between the bases no mat- 
ter how the linear combination of these bases are taken. 
This overlap is the origin of the ferromagnetic coupling 
between the electrons in the flat band. 



B. Known polymers of five-membered rings 

A typical polymer of five-membered rings is polypyr- 
role. The atomic structure and calculated band struc- 
tures are shown in Figs.[3fa) and (b), respectively. The 
five-membered rings alternate their directions in this 
compound, so that the unit cell contains two rings. X 
point in Fig. 0{b) nas k — (^/a, 0,0), where a is the 
lattice constant of this unit cell as shown in Fig. Eta). 
Because of the doubled unit cell, most of the bands are 
folded at X. For directions (y, z) perpendicular to the 
chain, we have taken a cell size large enough to avoid 
inter-chain interactions. We can see that polypyrrole has, 
not surprisingly, no flat bands around the Fermi energy. 

We have also looked at other typical polymers of five- 
membered rings, including polythiophen (where an N-H 
block in polypyrrole is replaced with an S atom) and 
polytriazole (where all the C-H blocks in polypyrrole are 
replaced with N atoms). However, their band structures 
have no flat bands, either. 



C. Designing the polymer 

This is exactly where the designing comes in. Our 
strategy consists of two approaches. One is to replace 
H atom bound to N atom at the top of the ring with 
various kinds of bases. According to the tight-binding 
calculation described in the previous subsection, the on- 
site energy (e\) of the top of the ring should be higher 
than that of C atoms to have a flat band. Because the 
on-site energy of N is expected to be lower than that of 
C, the first attempt should be the replacement of the N 
atom with other elements with higher on-site energies. 
However, that should be difficult if one wants to retain 
the existence of the 7r-orbital and single electron occu- 
pation on it. Therefore, we have opted for controlling 
the on-site energy by replacing the H atom with other 
functional groups, by leaving the N atom intact. 

Alternatively, we can replace C-H blocks on the bot- 
tom edge of the pentagons (3, 4, 8, 9 in Fig. Ufa)) w ith N 
atoms, where the lowering of the on-site energies of the 
bottom of the ring may effectively realize the required 
condition. In other words, we use polypyrrole and poly- 
triazole as the platform to modify either the top or the 
bottom of the ring. Since it is not obvious which ap- 
proach should be better, we consider both platforms for 
all the substituents considered here. 

Substituents we have tested are sodium (Na), potas- 
sium (K), chlorine (CI), fluorine (F), cyanogen (CN), 
nitro (NO2), sulfate (SO4), carboxyl (COOH), amino 
(NH 2 ), methyl (CH 3 ), and hydroxyl (OH). For all these 
candidates, we have performed first-principles (GGA- 
DFT) optimization of atomic structures, and obtained 
their band structures. Most of them have no flat bands 
around Ep, but some of them have turned out to have 
flat bands. For these we have performed further (GGA- 
SDFT) first-principles optimizations under the doping 
condition to make the flat band half-filled. 

First, we have examined single atoms (alkali metal or 
halogen) as substituents. This is motivated by the fact 
that small substituents may be better than larger ones, 
because a bulky group at the top of the ring may intro- 
duce a coupling (transfer) with the bottom of the adja- 
cent ring, which would destroy the basic assumption that 
the adjacent rings are connected by one bond to satisfy 
Mielke-Tasaki's condition. 

In order to raise the on-site energy of the top of the 
ring, we have first considered an alkali-metal atom, Na, 
because of its low electron affinity. Figure EJa) shows the 
atomic structure of poly-sodium-triazole, where hydro- 
gen atoms in polytriazole are replaced by sodium atoms. 
Figure |3Jb) shows the band structure, where one can see 
that a nearly flat band exists above the Fermi energy. 
Absence of any unnecessary bands between the flat band 
and Ep is good sign, because we can dope carriers into 
the flat band without any complications. 

We have also examined the Na-substituted polypyr- 
role, and obtained a similar band structure. However, in 
this case the flat band above Ep is much closer to a dis- 
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FIG. 3: (a) The optimized atomic structure of poly-sodium- 
triazole, the band structure of poly-sodium-triazole when un- 
doped (spin unpolarized)(b), and doped (spin polarized) (c). 

persive band below Ep, which can introduce unwanted 
inter-band interactions. As for the alkali-metal we have 
also considered potassium, but the case of Na has turned 
out to be better. 

We move on to the doping. First, the flat bands always 
appear as a pair of bands folded at X point, except for 
the case of antiferromagnetic ordering. Therefore, we 
make the flat-band half-filled by doping two carriers per 
unit cell into the pair of flat-bands. Namely, one carrier 
is doped for each ring. In poly-sodium-triazole, we have 
considered the doping of two electrons per unit-cell. Here 
we realize the doping condition by increasing the number 
of electrons with a uniform positive background charge 
for charge neutrality. 

After a geometrical optimization of doped poly- 
sodium-triazole, we have calculated the spin-dependent 
(GGA-SDFT) band structure in Fig.^c). One can see 
that the number of occupied states differs between major- 
and minor-spins, which means that the system is fer- 
romagnetic. It should be noted that (the pair of) the 



flat band is not entirely spin-split. The difference of the 
number of spins is 0.65 per unit cell. The spin polariza- 
tion is smaller than the expected polarization (2.0). We 
also found that the total energy (per unit cell through- 
out this paper) of the ferromagnetic (F) state is lower 
by 12 meV than that of the paramagnetic (P) state. We 
could not find antiferromagnetic (AF) solution even when 
we started geometry optimization from an AF electronic 
state. So the F state is the most stable in this material, 
but is only slightly lower in energy than the P state. 

The reason why this material gives such an insufficient 
result in spite of the existence of the nearly flat band 
should be sought in the nature of the flat-band's wave 
function. We found that the half-filled flat band here 
is made of an orbital localized around the Na atom, so 
does not satisfy the local connectivity condition. This 
should be why the stability of the F state is very weak, 
while the finite spin-polarization in GGA-SDFT should 
be a narrow-band effect rather than a much more robust 
Mielke-Tasaki effect. Indeed, the difference between the 
weak ferromagnetism here and a robust ferromagnetism 
in poly-aminotriazole, shown later, is a good demonstra- 
tion of the importance of the local connectivity condition, 
hallmark of the flat band ferromagnetism. 

Now we digress a bit, and consider the contrary case of 
an atom having high electron affinity such as a halogen, 
e.g., fluorine atom, which might be heuristic. Figurc^Ja) 
shows the atomic structure of poly-fluoropyrrole, where 
the fluorine atoms are substituted for hydrogen atoms 
bound to nitrogen atoms in polypyrrole. Figure Hfb) is 
the band structure obtained by spin-unpolarized (GGA- 
DFT) calculation. We find a flat band which lies just 
in between two dispersive bands below Ep. Because of 
this the doping greater than one carrier per unit cell is 
required. 

So we consider making the flat-band half-filled by dop- 
ing four holes per unit-cell. Here, two holes are necessary 
to make the dispersive band just below the Fermi level 
empty, and two holes are necessary to make the folded 
flat bands half-filled. Figure 0fc) shows spin-dependent 
band structure after the geometry optimization under the 
doping condition. The system is seen to be ferromag- 
netic. The difference of the number of spins per unit cell 
is just 2.0. One problem of this material, however, is 
that doping twice as large is necessary. When the dop- 
ing is performed in this material, the flat band shifts 
upward (downward) for minority (majority) spin. Re- 
placing the fluorine atom with a chlorine atom results in 
a band structure similar to the case of fluorine. In fact 
an AF state, obtained by starting geometry optimization 
from an antiferromagnetic state, has a total energy 21 
meV lower than that of F state. 

We move on to functional groups as substituents to 
search further possibilities. We have considered an exam- 
ple, poly-hydroxytriazole, where H atom in polytriazole 
is replaced by OH with a low electron affinity Although 
a pair of slightly dispersive bands exist just below Ep, 
these bands, when doped with two holes per unit cell, 
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FIG. 4: (a) The optimized atomic structure of poly- 
fluoropyrrole, the band structure of poly-fiuoropyrrole when 
undoped (spin unpolarized)(b), and doped (spin polar- 
ized) (c). 



become mixed with dispersive bands underneath. As a 
result, the difference of the number of spins in the F state 
is only 0.39 per unit cell, and the paramagnetic state is 
more stable than the F state by 3 meV. 

Figure [5fa) is the atomic structure of poly- 
aminopyrrole, where H atom bound to N atom in 
polypyrrole is replaced by NH2, whose electron affinity 
is low. In the band structure, Fig. Efb), one can find a 
pair of flat bands below Ep. The lower half of it is very 
flat while upper half is slightly dispersive. 

We have then doped two holes per unit cell. The 
band structure after geometry optimization is shown in 
Fig.[SJc). The difference of the number of spins is 1.0 per 
unit cell. As one can see from Fig. \^c) , the Fermi level 
intersects the middle of the pair of flat bands. Although 
the lower part of the pair of flat bands is quite flat, the 
dispersion of the upper half part is too large. Here only 
lower half part seems to work as a flat-band. The total 
energy of paramagnetic state in this material is found to 



FIG. 5: (a) The optimized atomic structure of poly- 
aminopyrrole, the band structure of poly-aminopyrrole when 
undoped (spin unpolarized) (b) , or doped (spin polarized) (c). 

be 65 meV higher than that of ferromagnetic state, while 
an AF solution does not exist, so the F state is the most 
stable. 

All other substituents tested here except for the one in 
the next section have turned out to be inappropriate. In 
most of those materials, the band which we expect to be 
flat is either too dispersive or too deep. For the remain- 
ing ones, band structures are qualitatively similar to the 
examples described in this section. Poly-aminotriazole, 
discussed in the next section, also gives qualitatively sim- 
ilar results to poly-aminopyrrole shown in this section, 
but quantitatively much better. 

III. POLY-AMINOTRIAZOLE 

A. First principles calculation 

The best candidate for ferromagnetic organic poly- 
mer among the materials we have considered is poly- 



(c) 



£ -8 



r x 

majority spin minority spin 



LU 



Ef 
-10 



-12 



-14 





x r 



x 



FIG. 6: (a) The optimized atomic structure of poly- 
aminotriazole, the band structure when undoped (spin unpo- 
larized)(b), and doped (spin polarized) (c). The bands having 
7r (ex) character are shown by solid (dashed) lines, respectively. 
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FIG. 7: A band structure for a metastable antiferromagnetic 
state in poly-aminotriazole. 



the number of spins is the desired 2.0 per unit cell. 

In this material, we found a metastable antiferromag- 
netic state. The band structure of the AF state is shown 
in Fig. □ This is calculated with GGA-SDFT, but the 
resulting band structures for up- and down-spins are the 
same, while the wave functions for up- and down-spins 
are arranged alternately on the chain of rings. This 
means that even though we have avoided intra-ring AF 
instability by choosing five membered rings, inter-ring 
AF can exist at least as a metastable state. The upper 
and lower parts of the flat band folded at X point is sep- 
arated with Ep in between. The AF state is found to 
be 52 meV higher in energy than the F state. Paramag- 
netic state is even higher (by 384 meV) than the F state. 
Therefore F state is the most stable in this material. 



aminotriazole2i. Figure EJa) shows the atomic struc- 
ture of poly-aminotriazole, where H atom in polytria- 
zole is replaced with NH2. The difference between poly- 
aminopyrrole shown in the last section and this material 
is whether the platform is polypyrrole or polytriazole, i.e., 
the C-H blocks in the bottom of the ring is replaced with 
N atoms. The band structure of the poly-aminotriazole, 
shown in Fig. E{b), is similar to that of Fig. EJb), but 
the dispersion of the upper half of the pair of flat bands 
is smaller and the separation from dispersive bands be- 
low is greater than that of Fig. Elb) . These features are 
desirable for the flat-band ferromagnetism. 

After the optimization of the atomic structure under 
the doping condition of two holes per unit cell, the spin 
dependent band structure was calculated and shown in 
Fig. [Sf c) . The result hits on the ideal situation — The 
pair of flat bands is made half-filled, and as a result, fully 
occupied by majority spins while totally unoccupied by 
minority spins. Reflecting this situation the difference of 



Let us discuss how the flat-band ferromagnetism is 
achieved in this material by combining the wave func- 
tions of the flat band obtained by first principles calcu- 
lation with the calculation in the tight-binding model. 
Figure IS] shows the majority-spin's wave functions in the 
pair of flat bands just below the Fermi level in Fig. EJc). 
While the phases of the wave functions can be taken ar- 
bitrarily, we choose for clarity the wave functions of the 
lower and upper bands at T point (Figs. EJa) and (f)) 
to be real and pure imaginary, respectively, and those at 
intermediate fc-points to be continuous. If we first look 
at the real part, typically for Fig.[HJa), the wave function 
for each monomer (aminotriazole) has two nodal lines 
which divide the amplitude into three parts, correspond- 
ing to NH2, C-N-C, and N-N blocks. For the imaginary 
part, typically for Fig. IBJf), the wave function has one 
nodal line, which divides the amplitude into two, con- 
sisting of two C-N blocks. The two wave functions (real 
and imaginary), which are orthogonal, are mixed at the 
intermediate fc-points. 
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FIG. 8: Wave functions of the polarized lower band (a, b, c) and upper band (d, e, f) of poly-aminotriazole. (a, f) for k x — 
(r), (b, e) for k x — ir/2a, and (c, d) for k x = ir/a (X). Inset at the left side indicates the band and fc-point of each figure. 




FIG. 9: (a) Band structure in the tight-binding model, (b) 
Tight-binding parameters used for the fitting, (c) The eigen- 
function which satisfies the connectivity condition. 



B. Tight-binding model 

In order to capture the essence of the flat band, here 
we introduce a tight-binding model which represents 7r- 
orbital's network in this system. Figure Ela) shows its 
band structure. The parameters of the model displayed 
in FigOJb) are fitted so as to reproduce Fig. HJb). Ob- 
tained values are: tcN — tec — U — 2.5 eV, iNN = 3 
eV, £o = —1-4 eV, and E\ = —0.5 cV. One can see that 
the tight-binding band (Fig.^Ja)) excellently agrees with 
the 7r-bands obtained by first-principles calculation (solid 
lines in Fig. EJb)). 

Figure |jjc) shows the tight-binding wave function of 
the eigenstate corresponding to that in Fig.^a). In this 
material, one site is added to the top of the ring which 
represents a nitrogen atom of NH2 . The open and closed 
circles indicate the sign of the wave function, respectively, 
and their sizes amplitudes. This wave function satisfies 
the local connectivity condition. 

Here it should be noted that the substitution of NH2 
for H not only raises the on-site energy of the top of the 
ring but also introduces an extra 7r-orbital of the N in 
NH2. The existence of the extra site completely changes 
the connectivity condition, so that we have no longer to 
satisfy e± — t > 0, which has turned out in the present 
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FIG. 10: fc-dependence of the wave functions in the tight- 
binding model. 



study to be rather difficult condition to realize. 

Figure shows /c-dependence of the tight-binding 
wave functions of the flat band. Here fc-space is extended 
to the second Brillouin zone, so the left and right halves 
of the figure represent the wave functions of the lower and 
upper parts of the pair of flat bands, respectively. Phase 
factors are chosen again for clarity as in Fig. [3] The la- 
bel "Re(2,5)" indicates that the line represents the real 
part of the wave function at the site 2 and 5, where the 
numbering of the sites is shown in the insets. On the T- 
points at the left- and right-edge of the figure, the wave 
function becomes real and pure imaginary, respectively. 
These two wave functions are displayed in the inset of 
Fig. ^3 where it should be noted that we display a part 
of the periodically extending Bloch functions, in contrast 
to Fig. c) which shows the entire eigenfunction. 

One can see that the symmetry and the character of 
the wave functions shown in Fig. 1101 is similar to the 
first-principles wave function (Fig. |HJ) . These wave func- 
tions can be constructed by a linear combination of the 
eigenfunctions shown in Fig. [IJc) . Figure ^2 illustrates 
how the linear combination of the localized eigenfunc- 
tions provides the Bloch function. When the two neigh- 
boring eigenfunctions, A and B, are added, the sum 
at the ring on which they overlap (enclosed by dot- 
ted lines) forms the shape of the imaginary-part wave 
function as shown in the inset of FigEH When sub- 
tracted, it forms the shape of the real-part wave func- 
tion. The entire imaginary- and real-part wave functions 
are constructed by the summation of all localized eigen- 
functions (A+B+C+D+- • •) and the combination of the 
summation and subtraction (A— B+C— D+- • •), respec- 
tively. The real-part wave function and the imaginary- 
part wave function are orthogonal to each other. The 
gradual mixing of these two orthogonal wave functions 
as k is changed is the very character of the Mielke-Tasaki 
flat band, which should never occur on an ordinary flat 
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FIG. 11: Schematical linear combinations of the neighboring 
eigenfunctions. 



band. 



C. Comparison with the Hubbard model 

Having constructed the tight-binding modeling we pro- 
ceed to the question of whether the ground state is spin- 
polarized in the presence of the Hubbard interaction, 
T-tu = UJ2i n i^ n il- I n Fig. El we show the phase dia- 
gram against U and S\ obtained with an exact diagonal- 
ization calculation for a 12-site (2 rings) Hubbard model 
for to — 2.5 eV and various values of £nn = 3.0 — 4.0 
eV. As indicated in the inset, £o is chosen throughout to 
satisfy the condition for the flat band, 

£ — £\ = (1 — e)(tNN — ( £ — £ l) 2 ) — *NN, 

(e — £\ + <nn)/(1 — e) = — £nn/(£o — e) — £nn(e — £o)j 

where s is the eigenenergy of the flat band and tcN — 
tec = to (= 1 here) is assumed for simplicity. We can 
see that we have indeed a ferromagnetic phase unless the 
repulsion is too strong (i.e., U < U c with U c = 2 ~ 5 eV). 
Our preliminary quantum Monte Carlo calculation con- 
firms that two unit cell is sufficient to roughly determine 
the ferromagnetic phase boundary. Peculiar destruction 
of the ferromagnetism above U > U c is another hallmark 
of the Mielke-Tasaki magnetism^ 



9 




fNN=3.0 



Ferromagnetic 

-°1.5 -1.0 -0.5 

£1 (eV) 

FIG. 12: Stability of the ferromagnetic phase calculated with 
an exact diagonalization of the Hubbard model. 



IV. DEVIATION FROM THE HALF-FILLED 
FLAT BAND 

A. Comparison of the total energies 

In order to realize ferromagnetism in the materials 
we proposed here, carriers have to be doped into the 
flat band. We can consider several possible methods to 
achieve this such as field effect dopingSS,, photo-induced 
carrier injection 2 ^, as well as chemical doping. The first 
two are now attracting much attention in molecular elec- 
tronics, but they are very recent techniques and may need 
further developments. We can instead concentrate on 
the chemical doping, where the polymers are crystallized 
along with the accepters, and we found it feasible. De- 
tails are presented elsewhere^. Now, if the flat-band 
ferromagnetism can occur not only just at the half-filling 
but also away from that, this would make the chemical 
doping much easier. 

We have to know how robustly the ferromagnetism sur- 
vives when the filling of the flat band deviates from just 
half-filled. Although there is a proofii and a numeri- 
cal studj*22. for a Mielke-Tasaki lattice which show that 
slight deviations do not wash out the ferromagnetism, 
those conclusions may not be simply applied to other 
models. For the present model, calculations taking full 
account of the electronic correlation such as the exact di- 
agonalization of the Hubbard model would be unsuitable, 
because changing the number of electrons in small sys- 
tems corresponds to too big a change in the filling. We 
have instead performed an SDFT calculation, because it 
can handle the fractional number of electrons (per unit 



cell) if a uniform background charge is introduced. 

At this point we can note the following. In SDFT cal- 
culations, electronic correlations are not fully taken into 
account, so we can question whether the flat band fer- 
romagnetism may be described within the SDFT frame- 
work. This problem has been addressed in our previ- 
ous publication 27 - where SDFT and exact diagonaliza- 
tion of a Hubbard model have given qualitatively con- 
sistent results. This is not surprising, since one of the re- 
markable features of the flat-band ferromagnetism is that 
the magnetism exists over the whole range (infinitesimal 
< U < oo ) of the interaction as proved rigorously, so 
the weak-U case (for which the methods such as SDFT 
is meant to be applicable) crosses continuously over to 
strong- U case. As far as the on-site repulsion is not too 
strong, the comparison between the gain of the exchange- 
correlation energy by spin polarization (F state) and the 
gain of the band energy by forming AF order (or equiv- 
alently, a spin density wave (SDW) with the periodicity 
of two five-membered rings) accompanied by a lattice re- 
laxation should predominantly determine which state is 
stable. These two energies can be accurately estimated 
by SDFT calculations. Therefore, we expect here that 
SDFT can serve our purpose, to see the band-filling de- 
pendence of the total energies of F, AF, and P states. 

Figure ED ( a ) shows band-filling dependence of the 
total energies of the F and AF states of the poly- 
aminotriazole measured from that of the P state. The 
position where the number of holes, rift, is 2.0 corre- 
sponds to the half-filling. The results for F, AF, and 
P are obtained by fully relaxing the atomic geometry for 
each electronic state. It can be seen that F state is stable 
in the range 1.1 < rift < 2.0, while P state is stable for 
rift < 1.1. In the region where F is stable, AF is found 
to be metastable. In the P-region, rih < 1.1, no AF state 
is found. This result shows that the ferromagnetism in 
this material is robust, even down to the quarter filling, 
against the deviation from the half-filling. The result 
that quarter-filling is acceptable greatly relaxes the diffi- 
culty in doping. 



B. Correlation with the atomic configuration 

Since the polymer proposed here has functional groups, 
we should also examine how the three-dimensional 
atomic configuration is correlated with the band-filling 
dependence. Namely, all the atoms do not lie on a plane: 
in poly-aminotriazole two hydrogen atoms in NH2 in 
the undoped polymer lie above the plane of the rings, 
while all the atoms are co-planar in the half-filled poly- 
aminotriazole. We can more closely look at the torsion 
angle, tp, defined as the angle with which C-N-N-H atoms 
are connected. Figure H"3T b) shows the torsion angles in 
the P and F states, respectively. From Figs. E2 a ) an d 
(b), one can see that the breakdown of the ferromag- 
netism with the decrease in the number of holes is corre- 
lated with the change in the atomic configuration. 
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FIG. 13: Doping dependence of (a) the total energy (relative 
to the paramagnetic state) of F(solid line), AF(dashed line), 
and F(fixed structure, dotted line, see text) states measured 
from that of P state, and (b) torsion angle of C-N-N-H in 
F (diamonds) and P (reversed triangles) states. Stable states 
are connected by a solid line, while metastable states by a 
dotted line. Insets show front, side, and top views of the 
atomic structure of the aminotriazole in P and F states. 



to the N atom and two H atoms are formed in the three 
directions of sp 3 , leaving one lone pair, NH2 has a three- 
dimensional structure. This time, the network consists 
of five 7r-orbitals of the rings and a sp 3 hybridized or- 
bital of the lone pair. This change is disadvantageous 
for the scheme of the flat-band formation (Fig. [5J due 
to the balance of the six it orbitals. This is considered 
why F state always takes a flat geometry while P state a 
non-flat geometry. 

To see what part of the energy affects the change in the 
stable state when rih is decreased, we have calculated the 
energy difference between the F and P states by fixing 
the atomic structure to the optimum one for the F state 
at rih = 2.0. The result is plotted as a dotted line in 
Fig. I13f a). This shows that the value is negative over 
the whole range plotted here, which indicates that the 
electronic part of the total energy always favors F state 
regardless of the filling of the flat band. Therefore, the 
change in the stable state at rih = 1-1 can be ascribed to 
the change in the structure, especially of NH 2 . If we can 
find a material in which there is no such change in the 
atomic structure but the parameter values are similar, 
we can expect an occurrence of the ferromagnetism even 
for a low-doping region less than quarter-filling. 

In summary, the present work describes an example of 
materials design from first principles for the purpose of 
realizing flat-band ferromagnetism in organic polymers. 
We have found a promising material, polyaminotriazole, 
for which why it is the best and the robustness of the 
ferromagnetism has been discussed. 
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